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ABSTRACT: Recently, Schmidt and Schnack (Physica A260 479 (1998)), following earlier 
references, reiterate that the specific heat of A non-interacting bosons in a one-dimensional harmonic 
well equals that of A non-interacting fermions in the same potential. We show that this peculiar 
relationship between heat capacities results from a more dramatic equivalence between bose and 
fermi systems. Namely, we prove that the excitation spectrums of such bose and fermi systems are 
spectrally equivalent. Two complementary proofs of this equivalence are provided, one based on an 
\ analysis of the dynamical symmetry group of the A-body system, the other using a combinatoric 
analysis. 
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I. Introduction: With the advent of dilute atomic BECt 1 ' 2 ' 3 ] and, recently, nearly degenerate 
dilute atomic fermi gas' 4 ' 5 ], there is renewed interest in understanding aspects of quantum many- 
body theory in inhomogeneous (in particular, harmonically trapped) systems. Since trapped, cooled 
atoms have properties that are, in principle, controllable to a degree unavailable in other systems 
(e.g., clusters and nuclei), they present new opportunities to study quantum mechanics and many- 
body theory. 

Although the ultra-cold dilute atomic gas systems are large compared to the coherence lengths, 
they are not homogeneous, due to the fact they are generally trapped in a (nearly) harmonic po- 
tential. In many of these systems the interparticle forces are significant. In this note, however, we 
ignore the interactions between atoms, with the motivation being to understand better the thermo- 
dynamic properties of N trapped non-interacting bosons and fermions. Recent work! 6 ' 7 ] describes 
strange relations between the equilibrium thermodynamics of these two systems. It was shown, for 
example, that the heat capacity (as a function of temperature) of N noninteracting bosons in a 
one-dimensional harmonic potential is the same as that of N noninteracting fermions in an identical 
potential. The respective partition functions for these systems are likewise closely related (see Ref. 

[7])- i „ ' 

These " coincidences" provide hints that a deeper underlying connection exists between bose and 
fermi gases in a harmonic well. In particular, the heat capacity and partition functions as functions 
of the inverse temperature (3 can be thought of as an "imaginary time" continuation of a fourier 
transform of the spectrum. The fact that the heat capacities are the same for all temperatures 
suggests that there should be a state-for-state, level-for-level correspondence between these non- 
interacting many-body bosonic and fermionic systems. We show that this is indeed the case, and 
below describe two independent proofs of the spectral equivalence of excitations in these systems. 
The first of these arguments relies on the dynamical symmetry group properties of these systems, 
while the other is based on combinatoric methods. 

II. Dynamical symmetry group approach: Consider N non-interacting bosons or fermions in 
a harmonic well. In either case, the energy differences between energy levels is an integer multiple 
of %u. Consequently, to analyze the spectral properties of each system it suffices to study the 
multiplicity per (total energy) level. Below we will show that the entire spectrum of both the 
bosonic and fermionic systems are isomorphic up to an overall energy shift. 

Classically, a system of iV non-interacting particles in a 1-d harmonic potential is identical to 
that of a single particle in an iV-dimensional isotropic harmonic potential. The system thus has an 
obvious spatial O(N) symmetry we call "angular momentum". However, it is apparent with more 
introspection that the system possesses a much larger dynamical symmetry group. Orbits in the 
iV-dimensional isotropic harmonic potential do not precess. In analogy with the Kepler problem, we 
say that there is a conserved Runge-Lenz vector (which may be thought of as the axis of the orbit 
in configuration space), and we thus expect the symmetry group to be enlarged. 

Since we will be interested in the quantization of the system, we describe the dynamical sym- 
metry enlargement through the operators of the associated quantum theory. To simplify notation, 
take hu = 1 throughout. Label the raising and lowering operators for the bosonic theory a] , with 
i = 1,...,N. The canonical commutation relations (for the bosonic case) are [cxi,at] = 5ij. The 

many-body hamiltonian operator of this noninteracting system is H = ^ aj a « + e, where e is an 
overall constant. 

We call the space of the eigenvalues of the ajcxj the state space. Equivalently, the state space 
is the integer lattice in the (+, ...,+) quadrant of iV-dimensional Euclidean space. State space is 
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not the Fock space, but is a useful auxiliary space from which we will construct the Fock space, 
and so we discuss it's properties. Let be orthonormal unit basis vectors in this Euclidean space 
associated with the eigenvalues of a\cii. We name several distinguished vectors in this space, namely 
the level vector k = and the root vectors li = ej — e^+i for i = 1, . . . , N — 1. We also define a 

spanning set of weight vectors via (r^, / j) = dij together with (r^, k) = 0. 

Note that the operators associated with l i: namely aja^ — a| +1 aj + i, are independent, and com- 
mute with each other (being all diagonal) and with the hamiltonian. The hamiltonian corresponds 
to the level vector. Furthermore, to each pair of particles / ^ j, there is an associated su(2) subal- 
gebra generated by {a\ai — cdaj, a\a,j + aja;, i{a\a,j — aja/)}. The application of the second or third 
operators in the above sit (2) subalgebra "shift" the first operator's eigenvalue by a combination of 
root vectors. Finally, note that the matrix of inner products = (hjj) of the root vectors is 
exactly the cartan matrix of su(N). Thus, we have identified the dynamical symmetry group of this 
system generated by the (trace-free part of the) products of a\dj to be su(N). 

We now construct the Fock space for both fermions and bosons from the state space by realizing 
the respective anti-symmetrizations and symmetrizations of the multi-particle Fock states as linear 
combinations of states in the state space that lie on the same Weyl group orbit. We make this cor- 
respondence precise with the following observations. Each state in the state space can be thought 
of as a particular product of single particle states, its coordinates (the components of the U are 
integers) are simply the harmonic oscillator level of each particle. Constructing the multi-particle 
state associated with that product of single particle states consists of combining all the states from 
single-particle label permutations. The permutation group Sn is generated by primitive transposi- 
tions (. . . , rii, rii + i, . . .) — > (. . . , rii+i, n^, . . .). Each of these primitive transpositions acts as a Weyl 
reflection (acting on all the roots) about the hyperplane perpendicular to the root h. 

Thus, the Weyl group, W of the symmetry algebra su(N) is exactly the group of permutations 
of the single particle states that make up the many-body state. Each element of the Weyl group 
preserves the level k. We specify a many body state through an assignment of a highest weight 
vector r and level s (a natural number) for which r + is a vector in the (+, ...,+) quadrant 
(boundaries included) . Explicitly, in terms of the vectors in the state space, the bosonic many-body 
Fock space has the basis ^ b r °l on 



where the sgn(a) is 1 if a is an even permutation and -1 if it is an odd permutation. Note that 
according to this definition, only r vectors from the interior of the Weyl chamber are associated with 
a fermionic many-body state. 

Succinctly stated, the multi-particle permutation symmetry of quantum mechanics maps the 
single particle states of state space into the highest weight space of the symmetry algebra su(N). 
For bosons, the map covers the entire Weyl chamber (including the lattice points in the bounding 
hyperplanes) at each level. For fermions, the map covers only the interior lattice points of the Weyl 
chamber. Additionally, due to the constraint that r + ^ is in the (+, ...,+) quadrant, at each level 
there are of course only a finite number of highest weight candidates. 




whereas the basis of the fermionic many-body Fock space is 
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The vector p = \ J2 a >o a (half the sum of positive roots) translates the vacuum of the bosonic 
Fock space to that of the fermionic Fock space at each level. Note also that p is thus orthogonal to 
the level vector k. It can be combined with the level vector to constitute a one-to-one map between 
the spectrum of the N boson and N fermion systems. Note that translation by the vector vector 
T = (0, 1, 2, . . . , N — 1) is precisely that map, and that T = p + ^-^-k. Note further that T has 
level r • k = N(N — l)/2, which is precisely the ground state energy shift between the bosonic and 
fermionic system. Geometrically, V is the smallest lattice vector that translates the lattice points 
in the bounding hyperplanes entirely into (the subset of) the interior of the Weyl chambers at each 
level. 



Ill Combinatoric approach: The spectral equivalence of one-dimensional, noninteracting 
harmonically-trapped bosonic and fermionic gases can also be understood through a straightfor- 
ward combinatoric argument. 

In a system of N noninteracting particles (bosons or fermions) in a harmonic well let the 
energy level of the zth particle be specified by the integer e;, with E = J2iLi e i the total energy 
of the system. (Note: In writing the energy as an integer, we are, as before, setting Tiuj = 1, 
and for notational convenience are ignoring the constant 1/2 associated with the single-particle 
ground state energy.) Clearly there are many different micro configurations possessing the same 
total energy E; we let Gn(E) denote the multiplicity of states with fixed energy E. We will show 
that the multiplicity functions for bosons and fermions are equivalent. More precisely, we show that 
G b £ son {E) = G f x rrnion (E + N(N - l)/2), indicating that the multiplicities for the bose and fermi 
cases are identical provided each is measured relative to its respective ground state energy (i.e., for 
bosons and N(N — l)/2 for fermions). This equivalence is sufficient to prove the spectral equivalence 
of the excitation spectrum. 

We begin with the bose case. We imagine ordering the iV particles from lowest energy to highest 
(ei, e2, . . • , ejsr). The energy of the lowest-energy particle (ei) can range from zero up to a maximum 
value of [E/N], where the brackets [ ] denote the integer part of the expression enclosed. (It is 
readily seen that if the energy of the lowest-energy particle were to exceed this maximum value, then 
the sum of the energies of the N individual particles would exceed the total specified energy E of 
the system.) 

For a fixed e\, the remaining energy E — e\ must be divided up among N — 1 particles. So the 
possible values of e2, which represents the lowest energy among the remaining (N — 1) particles, can 

range from ei to ^Erf • (As before, it is clear that if e2 went outside this range, then the sum of 
the energies of the N — 1 particles would exceed the prescribed value E — e\.) 
Proceeding in this fashion, we see that 

rB-ei 1 r- E - e l- e 2----- e JV-2l 

[E/N] [tv=t] I 5 J 

c b r n {E)= E E ••• E L 

ei=0 e 2 =ei ejv_i=ejv_2 



A similar argument is used to construct the multiplicity function for the fermionic case. The 
fundamental distinction stems from the additional constraint that two fermions cannot occupy the 
same energy orbital, which in turn modifies the lower and upper bounds in the above summations, 
as we now describe. Consider first the lower bounds. From the exclusion principle, it immediately 
follows that the lower (fermionic) bounds must take the form = e^-i + 1. The upper limits are 
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found by noting that for a system of N fermions with total energy E, the energy of the lowest-energy 
fermion cannot exceed 



N(N-l) 
2 

N 



as a straightforward calculation reveals. Consequently, we find 



s-ifermion 
N 



(E) 



N(N-l) 
2 

AT 

E 

ei=0 



(JV-l)(iV-2) 



E 

e 2 =ei + l 



(2)(1) 
'N-2- a 



E 

ejv_i=ejv-2 + l 



Expressed in this manner, the equivalence of G b £ son (E) and G f N ermion (E + N(N - l)/2) is now 
revealed through the following key coordinate transformation: In the fermionic summations above, 
introduce new coordinates = — i + 1. We claim that this will transform the fermionic sum 
into the corresponding bose sum. (Note: in the context of the preceding analysis, this coordinate 
change serves to relate the interior lattice points of the Weyl chamber (fermionic case) to the entire 
Weyl chamber (bose case), that is, it is simply the translation by the vector V) To see that this 
transformation achieves the desired result, first observe that under this transformation, the lower 
bounds in the fermionic summations (e$+i = e$ + 1) become (e^+i = e^), just as in the bose case. 
Meanwhile, it is not difficult to verify that the upper limits in the fermionic summations 



ej+i = 



E - e 1 - e 2 



(N-i)(N-i-l) 
2 



N -i 



now take the form 



ei+i 



E 



N(N-l) 
2 



ei - e 2 



N -i 



which, again is the same as for the bosonic case (once we shift by the fermionic ground state energy 
E -> E + N(N- l)/2). 

This equivalence between the bosonic and fermionic multiplicity functions proves that the ex- 
citation spectrum of one-dimensional harmonically trapped iV non-interacting bosons is identical to 
that of iV non-interacting fermions. 



IV. Remarks and Conclusion: Although the excitation spectra of the fermi and bose systems are 
identical, these systems are not related by an obvious supersymmetry. There may, however, exist a 
connection associated with the fermionic representation of affine lie algebra characters as described 
in Refs. [9,10,11]. Lastly, we observe that the recent work of Schmidt and Schnack 6 ' 7 indicates that 
the specific heats of similar bose and fermi systems in higher spatial dimensions (specifically, odd 
dimensions) might also be equivalent, just as for the one- dimensional case considered here. However, 
preliminary work suggests that spectral equivalence does not persist in higher (odd) dimensions. 
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